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The intensity statistics of signals in the presence of Gaussian noise is obtained by studying the 
model of a random signal plus a random phasor sum. The additive Gaussian noise is shown to 
result in a Bessel transform of the probability density of the signal intensity. The transformation of 
the intensity statistics can generally be applied to a mixture of independent random signals, one of 
which being a complex-valued Gaussian random process. It is used to retrieve intensity statistics of 
microwave pulsed transmission from Gaussian noise at long time delays. 
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I. INTRODUCTION 

The signal available at the output of a radio measure- 
ment system is never an entirely accurate indicator of the 
quantity to be measured. The accuracy of the measure- 
ment depends on the amount of noise compared to the 
level of signal, or signal-to-noise ratio. Signal noise can 
result from a variety of causes, both man-made and nat- 
ural. In most cases, however, it is natural additive noise 
which is the limiting factor for signal detection. This 
may include noise radiation from the sky picked up at 
the antenna, Johnson-Nyquist noise and shot noise gen- 
erated in the circuitry of the receiver, etc. Usually the 
natural additive noise can be represented mathematically 
as a Gaussian random process, hence the term Gaussian 
noise. The problems of detecting signals in Gaussian 
noise and of estimating parameters of signals in Gaus- 
sian noise have been studied by means of detection the- 
ory. There exists an extensive literature on this subject 
(see, for example, [l| for review). 

Our interest in the statistics of random signals in Gaus- 
sian noise arose from measurements of electromagnetic 
waves transmitted through random media. Wave trans- 
port in the presence of disorder can be characterized by 
the degree of nonlocal intensity correlation, which re- 
flects the closeness to the Anderson localization transi- 
tion (see, for example, and Q)- The presence of long- 
range correlation of intensity within a sample leads to en- 
hanced fluctuations of total transmission over the value 
predicted if correlation of intensity were short range, as 
is correlation of the field. The occurrence of enhanced 
transmission fluctuations can be seen in an ensemble of 
quasi-lD random samples, in which the sample length is 
much greater than the diameter of its cross section [4j] . In 
pulsed transmission measurements, the variance of trans- 
mission fluctuations normalized to the ensemble-averaged 
transmission increases with time delay from an exciting 
pulse Q while the decay rate of the average intensity 
within the sample decreases @ , reflecting two related ef- 
fects: the increasing impact of localization and the grow- 
ing weight of long-lived electromagnetic quasi-modes. At 
long times, however, the decaying transmitted intensity 
becomes comparable to Gaussian noise then affecting the 



intensity statistics so that the measured variance of total 
transmission is no longer the localization parameter. To 
study the dynamics of transport at long time delays, the 
impact of Gaussian noise must therefore be determined. 

To solve this problem, we consider the model of a ran- 
dom signal plus noise in which noise is represented by a 
random phasor sum with circular Gaussian statistics. We 
find that the addition of complex-valued Gaussian noise 
results in a Bessel transform of the probability density of 
the signal intensity. Depending on the physical problem 
under consideration, the solution can be used to find the 
intensity statistics of signal in the presence of Gaussian 
noise, to retrieve the intensity statistics of signal from 
Gaussian noise given the noise intensity, or to determine 
the noise intensity given the intensity distribution of the 
signal. We shall use it to determine the probability den- 
sity of pulsed transmitted intensity in Gaussian noise at 
long time delays. More generally, the solution can be 
applied to a mixture of independent random signals, one 
of which being a complex-valued Gaussian random pro- 
cess. An important example is the stationary field of 
a disordered cavity coupled to the environment, which 
can be represented by superposition of a standing wave 
(an eigenstate) and a traveling wave associated with the 
energy leaking out of the system |Tl'|, the former play- 
ing the role of 'signal' and the latter playing the role of 
'Gaussian noise'. 



II. STATISTICAL MODEL OF A RANDOM 
SIGNAL PLUS GAUSSIAN NOISE 



We consider the model of a random signal plus Gaus- 
sian noise, E — + En, m which Gaussian noise, E^, is 
represented by a random phasor sum with circular Gaus- 
sian statistics ['31 • The probability density of the real 
and imaginary parts of E^, and in, respectively, is a 
Gaussian with width a, 
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where (...) represents the average over an ensemble of re- 
ahzations. Because Es and are statisticahy indepen- 
dent, the probabihty density of the resultant field E is 
simply the convolution, P{E) = P{Es)*P{En). Finding 
statistics of the intensity of the resultant field, / = 
however, is more involved. 

The Gaussian statistics of and results into the 
exponential probability density of the noise intensity. 
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where (/„) — 2(7^ = D is the average noise intensity that 
we shall denote by D. The moments of the noise intensity 
follow from (3) as (J") = nl D". The signal intensity. Is, 
can be written as Is = |-Esp = + i^, where Ts and is 
are the real and imaginary parts of Eg, respectively. The 
moments of Is are given by 
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where P{Is) is the probability density of the signal inten- 
sity. 

To find the probability density of the intensity of the 
resultant field, we first calculate its moments (/"). Ex- 
pressing /" in terms of the real and imaginary parts of 
the signal and noise, taking the average and using (2), 
we arrive at 
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From (5), the average and the variance are (/) = {Is) +D 
and var(/) = var(/s) -|- 2{Is)D + D^, respectively. 

From the moments (/") the characteristic function and 
probability density of I can be obtained Here these 
are derived using the characteristic function M/, of the 
signal intensity Is, 

MiSp) ^ (eM-pQ) - E ^(^s">- (6) 

n=a 

The moments (/") can be deduced from (5) and written 
as 
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where i„ is the Laguerre polynomial of order n. Sub- 
stituting (7) into (6) and by making use of a generating 
function of Laguerre polynomials |8[ , we obtain 
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, \pD\ < 1. (8) 



By changing variables, pD = sD/{l + sD), Eq. (8) can 
be written as 
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where M/(s) is the characteristic function of the intensity 
/. P{I) is related to Mi{s) in the usual way, 

/ioc J 
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Substituting Mj{s) of (9) into (10), we obtain 
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Changing the order of integration and integrating over s, 
we arrive at 
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where Iq is a modified Bessel function of the first kind of 
zero order. Eq. (12) is the main result of the paper and 
represents the transformation of the probability density 
of the signal intensity in the presence of Gaussian noise 
with the average intensity D. Some examples are in or- 
der. For a constant signal with intensity A, for which 
P{Is) ^ S{Is - A), we obtain from (12) 
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which is in agreement with the result of the model of a 
constant phasor plus a random phasor sum 0, [13) Ell- 
in the case when Es is itself a random phasor sum, 
P{Is) = exp(-/s/(/s)) /(/g). Then, as may be ex- 
pected, Eq. (12) yields P(/) = exp(-//(/)) /(/), where 
(/) = {Is) + D. Finally, in the case of a disordered 
cavity in the absence of non-proportional damping , 
P{Is) = cxp(-/s/2(JsyWi/27r/s(/s), which is the Porter- 
Thomas distribution Jll'|. We then obtain from (12) the 
probability density, P{x — I /{I)), in the crossover from 
closed to open system. 
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where b = {2D{Is) + D'^)/{{Is) + D)"^ , in agreement with 
Ref. [nl. 

The equation (12) can be inverted to yield the proba- 
bility density of the signal intensity. 
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However, Eq. (15) is not particularly useful to retrieve 
P{Is) from Gaussian noise with the average intensity D, 
because it requires a continuation of P{I) to negative /, 
which is not available from the measurement. Instead, 
P{Is) can be found by solving the equation (12). In 
the next section, we determine the probability density 
of pulsed transmitted intensity in Gaussian noise at long 
time delays. 
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III. INTENSITY STATISTICS OF PULSED 
TRANSMISSION IN GAUSSIAN NOISE 



Here the resuhs of statistical model of the previous 
section are used to retrieve the intensity statistics of 
pulsed microwave transmission through random media 
from Gaussian noise at long delay times 0, Q . Accurate 
measurements of the time-resolved statistics of pulsed 
transmission are absolutely essential for systematic study 
of wave transport in the presence of disorder. Spectral 
measurements of the field transmission coefficient of mi- 
crowave radiation were made in an ensemble of random 
dielectric samples, as described in The response to 

a pulse with a Gaussian temporal envelope at carrier fre- 
quency vq is obtained by Fourier transforming the prod- 
uct of the field transmission spectrum and a Gaussian 
spectral function of width Sv. The field of the temporal 
response is squared to give the intensity I(t) for each sam- 
ple configuration. The average intensity {I{t)) is found 
by averaging over the ensemble of realizations. The re- 
sult is shown on a logarithmic scale as the black solid line 
in Fig. la. The noise in the transmitted field manifests 
itself as a constant background in Fig. la. The analysis 
of the probability density i) in the negative time be- 
fore the pulse, i.e., when E = En, shows that P{rn, in) is 
a circular Gaussian. The average pulsed transmitted in- 
tensity (/s(i)) is then (/s(i)) = {I{t))-D, where D is the 
constant background in Fig. la. Once this background 
is subtracted (blue dashed line), the dynamic range is 
significantly enhanced. 

The variance of normalized intensity, var[/(i)/(/(t))], 
is shown as the black solid line in Fig. lb. In the 
pulsed measurement the variance of normalized trans- 
mitted intensity is expected to increase with time de- 
lay from an exciting pulse [H] as increasingly more of 
energy within the medium is stored in long-lived local- 
ized modes. In Fig. lb, in contrast, the variance of 
normalized intensity falls at long times to a value of 
unity, reflecting the increasing impact of Gaussian noise. 
The variance of transmitted intensity follows from (7) 
as var[/s(i)] = var[/(t)] - 2{I{t))D + D^. The evolution 
with time of the normalized variance, var[/s(i)/(/s(t))], 
is shown as the blue dashed line in Fig. lb. 

As it follows from Fig. 1, the probability distribution 
of the transmitted intensity is increasingly affected by 
Gaussian noise at long times. The probability density 
P{I/{I)) for the time delays t = 338, 364, and 389 ns, 
indicated by vertical dashed lines in Fig. 1, are shown as 
the black solid curves in Fig. 2. At these time delays, the 
relative amount of noise is D/{I) = 0.30, 0.51, and 0.72, 
respectively. Also shown in Fig. 2 is the (exponential) dis- 
tribution of the noise intensity found in the negative time 
(t = —224 ns). Apart from the uppermost curve, each of 
the curves is displaced by a multiple of two decades for 
clarity of presentation. In order to find the probability 
density P(/s/(/s)), Eq. (12) is to be solved. However, the 
form of P{Is/{Is)) is already known 0,i,[ll,[ll, 
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FIG. 1; (color on line). Measurements of the time-resolved 
statistics of pulsed microwave transmission in an ensemble of 
random dielectric samples. The average (a) and the variance 
normalized to the average (b) of the measured intensity I{t) 
(black solid lines) and the transmitted intensity Is{t) (blue 
dashed lines) deduced from Eq. (7). 
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and K — (var[/s(t)/ (/s(i))] — l)/2, that is, the probability 
density of the normalized transmitted intensity is given 
in terms of a single parameter, its variance. The values 
of var[/s(i)/(/s(i))] corresponding to the three time de- 
lays are 2.48, 2.82, and 3.42, respectively, as found from 
Fig. lb. The respective probability densities P(/s/(/s)) 
are deduced from Eq. (16) and shown as the blue dashed 
curves in Fig. 2. Thus we found P(/s/(/s)) without solv- 
ing Eq. (12). To check the validity of the P(/s/(/s)), we 
use Eq. (12) to deduce P{I/{I)) and compare it to the 
measurement. The calculated P{I / {!)) are shown as the 
red dotted curves in Fig. 2 and are in excellent agreement 
with the measured densities. 



4 



>i 

'u5 
c 

CD 
T3 
>, 

CO 

o 



10" 
10-^ 
10-^ 
10-3 
10-^ 
10-= 
10-^ 
10-^ 
10-3 



10-' 
10 

10- 

10 



10 



12 



1 1 1 1 1 




1 1 1 1 1 1 1 1 1^ 


r\ 




t=338ns ^ 


- \- 
r \ 




t=364ns_ = 


^ \ 




i=ooyns 




\ t= -224ns 




- , , , , 1 


, , 1 , , ,',\ , , , , 


, , , , 1 , , , , - 



IV. CONCLUSIONS 



In conclusion, we have found the transformation of the 
intensity statistics of random signals in the presence of 
additive Gaussian noise. The transformation of the inten- 
sity probability density is given by the Bessel transform 
of Eq. (12). This can be solved to retrieve the intensity 
statistics of signal from Gaussian noise, given the average 
noise intensity. We used Eq. (12) to determine the inten- 
sity statistics of pulsed microwave transmission through 
random media from Gaussian noise at long delay times. 
More generally, the results of this work can be applied to 
a mixture of independent random signals, one of which 
is a complex-valued Gaussian random process. 
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FIG. 2: (color on line). Probability densities of the pulsed 
transmitted intensities at the time delays t = 338, 364, and 
389 ns, and in the negative time {t = —224 ns), indicated by 
the vertical dashed lines in Fig. 1. The black solid curves are 
the measured probability densities P(//(/)). The blue dashed 
curves are the probability densities P(/s/(7s)), deduced from 
Eq. (16) using the corresponding values of var[Js(t)/{/s(t))] 
found from Fig. lb. The red dotted curves shown through 
the data are P(I / {!)) deduced from Eq. (12). 
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